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Abstract. The collective behaviour of statistical systems close to critical points is characterized by an 
extremely slow dynamics which, in the thermodynamic limit, eventually prevents them from relaxing to an 
equilibrium state after a change in the thermodynamic control parameters. The non-equilibrium evolution 
following this change displays some of the features typically observed in glassy materials, such as ageing, 
and it can be monitored via dynamic susceptibilities and correlation functions of the order parameter, the 
scaling behaviour of which is characterized by universal exponents, scaling functions, and amplitude ratios. 
This universality allows one to calculate these quantities in suitable simplified models and field-theoretical 
methods are a natural and viable approach for this analysis. In addition, if a statistical system is spatially 
confined, universal Casimir-like forces acting on the confining surfaces emerge and they build up in time 
when the temperature of the system is tuned to its critical value. We review here some of the theoretical 
results that have been obtained in recent years for universal quantities, such as the fiuctuation-dissipation 
ratio, associated with the non-equilibrium critical dynamics, with particular focus on the Ising model with 
Glauber dynamics in the bulk. The non-equilibrium dynamics of the Casimir force acting in a film is 
discussed within the Gaussian model. 

PACS. 64.60.Ht Dynamic critical phenomena - 64. 60. an Finite-size systems 



! 1 introduction 

Critical points occur in the phase diagrams of a variety 
of microscopically different systems, ranging from mag- 
netic materials (ferromagnetic/paramagnetic transition), 
, to pure fluids (liquid/ vapour critical point), binary mix- 
• tures (mixing/demixing transition) and strongly interact- 
' ing matter (QCD). In spite of these differences at the mi- 
, croscopic scale, a unified picture of the collective critical 
' behaviour emerges in terms of the so-called order param- 
, eter if of the transition, the nature of which depends on 
I the specific system while its fluctuations determine the 
physical behaviour close to critical points. In particular, 
in the case of ferromagnetic materials, </3(x, t) can be iden- 
tifled with the coarse-grained magnetization density at 
point X and time t, whereas for a fluid A/B binary mix- 
ture with critical concentrations c"'*/cg'', ip is given al- 
ternatively, by CA — c^'' or c_B — eg'* where ca,b are the 
space-time dependent concentrations of the species A and 
B, respectively. The thermal fluctuations of the order pa- 
rameter field (fi are correlated in space and time across 
a typical correlation length ^ and a (linear) relaxation 
time t}i. Upon approaching the critical point, ^ and tji 
become much larger than the corresponding microscopic 
length and time scales ^micr and Tmicrj respectively, result- 
ing in a spatially collective and temporally slow behaviour 
which is characterized by a certain degree of universality. 



In fact, the physics at scales much larger than Tmicr and 
^micr becomes largely independent of the microscopic de- 
tails of the system, depending only on its gross features 
such as the range of the microscopic interactions, internal 
symmetries, spatial dimensionality d, conservation laws 
etc. which characterize the so-called universality class of 
the transition. Universality is primarily an experimental 
fact based on the evidence that, e.g., ^ ^ £,o\r\~'^ and 
tfl ^ rolr]^"^ for r = {T - Tc)/Tc 0, where Tc is the 
critical temperature, and tq specifically depend on the 
microscopic details of the system and therefore are non- 
universal, whereas v and z are universal exponents in the 
sense specified above |l|2j . Universality can be exploited in 
order to study the critical behaviour of a statistical sys- 
tem via suitable minimal models having the same gross 
feature as the actual system and therefore belonging to 
the same universality class. In most of the cases this is 
done in terms of effective theories for the order parame- 
ter field <^(x, t) and the resulting critical scaling properties 
can be analyzed by means of powerful field-theoretical and 
renormalization-group methods [1] . 

The slow and collective behaviour which emerges close 
to critical points affects relaxation phenomena and the 
way in which they actually occur. In Sect. [2] we discuss 
the case of bulk systems, with particular emphasis on the 
emerging non-equilibrium properties, whereas in Sect. [3] 
we focus on confined critical systems in film geometry. 
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2 Non-equilibrium relaxation in the bulk 

Altough the picture presented below is quite general, we 
shall refer specifically to the case of uniaxial ferromagnets 
belonging to the Ising universality class, with a critical 
temperature Tc and a space-time dependent order param- 
eter ip{'x.,t) (magnetization density). 

One of the simplest instance of relaxation phenomena 
can be observed after a sudden thermal perturbation. In 
fact, assume that a system, initially in equilibrium with a 
thermal bath at temperature Tq, is brought at time t — 
in contact with a thermal bath at temperature T ^ Tq. 
As a consequence of this perturbation the system starts 
evolving out of equilibrium in a way which is expected to 
depend on the specific initial condition (e.g., on Tq). Af- 
ter this transient regime, of typical duration tEQ{T), the 
system equilibrates with the thermal bath at temperature 
T, loses memory of the thermal quench and its dynamics 
becomes stationary and invariant (in the absence of exter- 
nal fields) under time reversal; fluctuations are therefore 
described by equilibrium dynamics. The duration t^ciiT) 
of the transient regime, i.e., of the relaxation, depends on 
the temperature T and is such that t^qiT > Tc) is finite 
whereas tEQ(T < Tc) is actually infinite (in the thermo- 
dynamic limit) and the non-equilibrium relaxation goes 
on forever. In particular, we focus here on the case of a 
quench to the critical point T ^ Tc, referring the reader 
to Ref. [3j for a discussion of the phase ordering dynamics 
for T < Tc. The fact that tsqiT = Tc) is not finite is 
actually due to the emergent slow and collective critical 
behaviour. In this context it is important to understand 
the general features of this non-equilibrium relaxation, ex- 
ploring, e.g., how it depends on the initial condition of the 
system. In particular - referring to Ising-like magnets - we 
will consider the cases depicted in fig. [T] of quenches from 
an initial state with vanishing (A) and non-vanishing (B) 
magnetization. In addition to the dependence on the ini- 
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Fig. 1. Quenches in the phase diagram of an Ising ferromagnet: 
T is the temperature of the bath, H the external magnetic field 
and (T, H) = {Tc, 0) is the critical point. The initial conditions 
A and B are characterized by a vanishing and non-vanishing 
magnetization, respectively. 



tial condition, it is also of interest to understand how the 
universality which characterizes critical phenomena man- 
ifests itself in the long-time relaxation and if a scaling 
behaviour actually emerges. Some of these issues were al- 
ready addressed more than 30 years ago, with particular 
focus on the dynamics of the order parameter [415] . Only 
quite recently, however, it has been realized [5] that dy- 
namic correlations and susceptibilities display very inter- 
esting non-equilibrium properties such as those observed 



in glassy systems, for example ageing (see, e.g., Ref. [7] 
for a review). In this respect, critical dynamics provides 
a valuable, non-trivial but relatively simple instance of 
dynamics with ageing, which can be characterized to an 
extent actually out of reach in more general cases of glassy 
dynamics. 



2.1 Dynamic observables 

The relaxation process in a system of large volume V can 
be monitored by looking at the global order parameter 
<P{t) = V~^^^^y(p{x,t), its average over the possible 
realization of the thermal noise M{t) = {^{t)) and its two- 
time connected correlation function C{t, s) = (<?(t)??(s))conn- 
In addition to C{t, s) it is useful to analyse the linear re- 
sponse R{t, s) = 6M{t) / 6H{s)\h=o of the global order pa- 
rameter M{t) to a magnetic perturbation H{s) applied at 
time s < t. With the aim of highlighting non-equilibrium 
behaviours, we consider the so-called fluctuation-dissipation 
ratio (originally introduced in the study of glasses [8l6j ) 



Xit,s] 



_ TR{t,s) 



dsC{t, s) ' 

and in particular its limit for well-separated times 



lim lim X{t, s) 

t — ^oo s — *oo 



(1) 



(2) 



Indeed, for t,s^ tEQ{T) the fluctuation-dissipation the- 
orem (see, e.g., Ref. [7]) implies X{t,s) = 1 and there- 
fore X°° = 1 whenever equilibrium is eventually attained, 
which is always the case for T > Tc, X°° ^ f, instead, 
signals that the system evolves out of equilibrium even at 
long times. In particular, for quenches from the disordered 
initial state {M{t < 0) = 0) to T < and H = (do- 
main coarsening), scaling arguments yield X°° — 0, lead- 
ing to the global picture in Fig. [2] ^ . The value of for 



FDT^ 1 - 
- 



equilibr. 



Fig. 2. X°° as a function of the final temperature T of the 
quench starting from a disordered initial state W . The value of 
X°° for T — Tc depends on the universality class of the system. 



T = Tc, on the other hand, cannot be fixed on the basis of 
general arguments but it depends on the specific system 
(actually on the universality class which it belongs to [S])- 
Note that, in contrast to the case of the non-equilibrium 
dynamics for T < Tc (phase ordering) in which thermal 
fluctuations do not affect the leading scaling behaviour [3] , 
the dynamics at T = Tc is actually driven by them and 
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by the emerging collective behaviour of the system. In re- 
cent years has attracted a lot of attention in connec- 
tion to the fact that, formally, it might be used to define 
an effective non-equilibrium temperature T^s for the sys- 
tem (see, e.g., Ref. and references therein). Indeed, at 
large times, the rhs of Eq. ([T]) equals 1 (as required by 
the equilibrium fluctuation-dissipation theorem) if one as- 
sumes that the temperature of the system is not given by 
the temperature T of the bath but by Tcs = T/ X°° . In the 
case of mean-field glassy models it has been shown that 
Teff is actually a bona fide thermodynamic temperature 
which controls, e.g., the direction of heat flows [7]. In this 
context, the non-equilibrium critical relaxation provides a 
simple instance of slow dynamics in which such a property 
can be tested beyond mean field [TIT. 

To give answers to the various questions posed so far, 
it is possible to fully exploit universality [911 Ij in order 
to characterize within a field-theoretical approach J12'll] 
the non-equilibrium critical relaxation and the associated 
ageing phenomena (see Ref. [13] for a pedagogical intro- 
duction). This approach allows a systematic analysis of 
several aspects of these relaxation phenomena and yields 
analytic predictions for scaling functions, exponents, and 
amplitude ratios which characterize the scaling behaviour 
of correlation and response functions at large times, within 
different universality classes and dynamics. In addition, 
it is possible to highlight dynamical crossovers in C(i, s), 
R{t, s) [14 15J and in the persistence properties oiM{t) [l6] . 
due to different initial conditions. Within this field-theoret- 
ical approach, for example, instead of studying on a lattice 
'L'^ the non-equilibrium dynamics of the Ising model with 
spin-flip Glauber dynamics (which captures the behaviour 
of some anisotropic ferromagnets and alloys), one studies 
the Landau-Ginzburg effective Hamiltonian H-lip] [T for 
the order parameter (p(x, t) (x e R'') with a suitable re- 
laxational dynamics, known as Model A [2j: 



dt<p{^,t) - -D 



Sip{x, t) 



C(x,t). 



(3) 



D is a kinetic coefficient and C(x, t) a zero-mean stochas- 
tic Gaussian noise with (C(x, i)(^(x', i')) = 2A;bTD5(x — 
x')5(t — t'). The results presented below refer to this case. 



2.2 Scaling behaviour 

The renormalization-group analysis of Model A predicts 
that R{t,s) and C{t,s) display the following scaling be- 
haviours after a quench to T = Tc |17j : 



R{t,s)^ARt^ (-) TR{s/t,t/to), 
s<t, { y{^e-i (4) 

where a = (2 — rj — z)/ z is given in terms of well-known 
universal equilibrium static (ry) and dynamic (z) expo- 
nents [I], is the non-equilibrium universal initial-slip 



exponent [17J, J-^r,c are universal scaling functions nor- 
malized such that .?^fl,c(OjO) = 1. Ar c are non-universal 
constants whereas to is a non-universal time scale set by 
the initial value of the magnetization Mq = M(t = 0) and 
displays a universal dependence on it 



BmM, 







(5) 



where the universal scaling exponent k > is given, in 
terms of static and dynamic equilibrium and non-equili- 
brium exponents, by k = 6+a+P/ (vz). The non-universal 
amplitude Bm can be fixed by suitable normalization con- 
ditions (see Refs. |17ll4j for details). 

Interestingly enough, Eqs. ([4]) and (O indicate that 
of the initial condition only the value Mq of the magne- 
tization really matters in determining the scaling prop- 
erties of the ensuing relaxation. A more detailed analysis 
shows that correlations in the initial state (as long as they 
are short ranged) do only contribute to corrections to the 
scaling behaviour [17] . According to Eq. (|4]) , the two-time 
quantities C{t, s) and R{t, s) are homogeneous functions 
of the three time scales t, s, and tg. In particular, when 
s < t ^to, which is always the case if Mq = 0, the scaling 
form of R (analogous one for C) becomes 



R{t<^to,s) ^ ARt'^it/sff'i^^s/t), 



(6) 



where fR\x) — Tr^XjO). Equation ^ clearly displays 
the scaling behaviour typical of ageing phenomena: As a 
function of the time t at which the effect of the mag- 
netic perturbation is measured, the relaxation time is set 
by the time s at which the perturbation was applied, 
which is also referred to as the age of the system, be- 
ing the time elapsed since the quench. Upon increasing 
s the "response" to the magnetic perturbation becomes 
increasingly slow. The scaling behaviour in Eq. ([6]) was 
already spelled out in Ref. [T7] even tough the connec- 
tion with ageing has been realized some years later [Hj. 
If Mq = the scaling in Eq. ^ (and analogous for C) 
is valid also at long times and allows one to express X°° 
as an amplitude ratio = Ar/[{1 - 6) Ac] [MTT] . This 
amplitude ratio turns out to be universal and therefore 
its value calculated within the field-theoretical approach 
can be compared with the corresponding results obtained 
on the basis of different models in the same universality 
class (e.g., Ising model with Glauber dynamics studied via 
Monte Carlo simulations). 

In the opposite limit of large times compared to to, 
i.e., to <C s < t the scaling form of R (analogous one for 
C) becomes jT4] : 



R{t,s > to) 



aRt^{t/s)-P'l'^^^-^f^^: 



\slt), (7) 



where fR°\x) ^ x^'^^^/^^'^'' !FR{x,y oo) and ur is de- 
termined such that /jj°°''(0) = 1. Equation ([7]) displays 
a scaling behaviour (with ageing) analogous to the one in 
Eq. ([S]), with the major difference that the non- equilibrium 
exponent 9 has been replaced by the combination —l35/{vz) 
of equilibrium static and dynamic exponents. Also in this 
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case X°° can be expressed as an amplitude ratio X°° — 
afl/{[l + (36/ {iyz)]ac} where ac is the non-universal con- 
stant in the scaling form of C which corresponds to aji. Ac- 
cording to Eqs. ([6]) and ([7]), as soon as the initial value Mq 
of the magnetization is non- vanishing (quench from point 
B in Fig. [1]), the asymptotic behaviour for large times 
(and therefore the value of X°°) changes compared to the 
case in which the initial value of the magnetization is zero 
(quench from point A in Fig.[T]). An analogous conclusion 
has been drawn for different universality classes |18|14|15] . 
The full scaling functions in Eq. ^ actually describe the 
crossover which occurs, for a finite and fixed value of the 
magnetization Mq (and therefore to) when both times s 
and t increases from the case s < i ^ to (described by 
Eq. ([6])) to the asymptotic case to ^ s < t (described by 
Eq. dZD). 

The field-theoretical approach yields not only the pre- 
dictions, previously discussed, for the scaling behaviour 
and the associated exponents but also analytic expres- 
sions for other relevant universal quantities such as the 
fluctuation-dissipation ratio and its asymptotic limit X°° . 
Here we focus on the specific case of the universality class 
of the Ising model with Glauber dynamics, although anal- 
ogous predictions are available for different universality 
classes and dynamics (see Ref. [H] for a review). In Fig. [3] 
we report the dependence of the asymptotic fluctuation- 
dissipation ratio (see Eqs. dU) and ^) on the space 
dimensionality d of the model, after a quench from an 
initial state with Mo = (a) and Mq 7^ (b) to the 
critical point. The solid lines correspond to different Fade 
approximants of the O(e^) (a) [Hj and 0(e) (b) [13] series 
which have been calculated within the field-theoretical e- 
expansion, where e = 4 — d. For d > 4 the mean-field pre- 
dictions X°°{Mo = 0) = 1/2 and X°°{Mo ^ 0) = 4/5 
become exact. The available Monte Carlo estimates in 
d = 2, 3,4 [9l20] for Mq = and in d = 2 for 7\fo ^ [H], 
reported in Fig. [3Ka) and (b), respectively, are in quite 
good quantitative agreement with the analytical predic- 
tions. 

In summary, field-theoretical methods are a viable ap- 
proach to investigate non-equilibrium and ageing phenom- 
ena during critical relaxation in the bulk. The long-time 
properties of the non-equilibrium relaxation turn out to 
depend on the initial condition (e.g., Tq) via the corre- 
sponding magnetization Mo and novel crossovers in the 
response, correlation functions and persistence properties 
of M{t) [1^ occur for finite Mq, as confirmed by Monte 
Carlo simulations. Within this approach one can also in- 
vestigate the dependence of Ar°° ~ and therefore of the 
effective temperature Tcs — T/X°° - on the observable 
used for its definition. Interestingly enough, one finds TU] 
that X°° (equivalently, Tes) is independent of it only if 
the fluctuations in the system are of Gaussian nature (as 
in mean-field or spherical models) . This observation might 
be of relevance to the case of more complex glassy systems. 

In addition to the case illustrated in this Section, a 
variety of different universality classes and dynamics, both 
in the bulk and close to surfaces of semi-infinite systems 
have been investigated in the literature [21|llj . 



3 Relaxation in confined geometries 

In Section [2] we discussed some aspects of critical relax- 
ation in bulk systems, with particular focus on the non- 
equilibrium character of such a dynamics and on quan- 
titative predictions for the Ising universality class with 
purely relaxational dynamics. In this Section we consider, 
instead, relaxation phenomena in confined geometries, fo- 
cussing on fluctuation-induced forces and on the way they 
build up in time after a thermal quench. Perhaps, the most 
widely known example of such forces is the Casimir effect 
in quantum electrodynamics (QED) [22j : less known is 
the analogous effect which occurs in statistical physics, 
discussed for the first time by Fisher and de Gennes |23j . 
Indeed, whenever one confines a medium in which fiuctu- 
ations (of different physical nature) take place, effective 
forces arise on the confining walls. The medium can be 
constitued by, e.g., ''^He, ^He/'^He mixtures, classical bi- 
nary mixtures, a Bose gas or even a magnetic material. 
At variance with QED, in statistical physics such effective 
forces have generically a finite range which is related to the 
typical correlation length ^ of the confined fiuctuations. 
When this correlation length becomes much larger than 
the typical microscopic length scale ^micr in the system, 
the Casimir-like force Fc acting on the confining walls 
takes a universal form 

^ - '-^WO, (8) 

where A is the large transverse area of the d— 1-dimensional 
walls, assumed to be parallel and separated by a distance 
L. is a universal scaling function and therefore it is ac- 
tually determined by the gross features of: (i) the system 
in the bulk and of the fiuctuating order parameter, as 
mentioned in the Introduction, (ii) the gross features of 
the surfaces, which - to some extent - are summarized 
by the boundary conditions they impose on the order 
parameter (p{x,t) (see, e.g., Ref. [24j for a review) and 
(iii) the geometry of the boundaries. In what follows we 
shall focus on the case of parallel and infinite confining 
walls, i.e., on the film geometry. Due to universality, any 
minimal model with the same gross features as the con- 
fined systems we are interested in captures its universal 
behaviour, including the Casimir force Fc- This model 
is usually specified in terms of an effective Hamiltonian 
Ti.[ip] which determines the equilibrium distribution func- 
tion -Pcq[</'] ^ 6xp{— 7i[(p]/(fcBT)} of^ the order parameter 
and therefore the effective free energy of the model. In 
turn, close to the bulk critical temperature Tc and for a sys- 
tem large enough (i.e., L » Anicr), J' decomposes as the 
sum of a bulk term proportional to the volume V = Ax L 
of the system, one proportional to its surface area A and a 
third one which yields the leading and universal finite-size 
correction we are interested in: 

T^^kBTcAx 

r 1 1 (9) 

i/bulk(r) + /surf (T) + j^O{L/0 + ■■■ , 

where higher-order terms in L^^^^^^^ x (i'micr/-^) have been 
neglected. As a consequence of this decomposition, the 
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Fig. 3. Asymptotic fluctuation-dissipation ratio (see Eqs. ([TJ and ([2])) as a function of the space dimensionality d, within 
the Ising universality class with purely relaxational dynamics after a quench from a disorder (a) and magnetized (b) initial 
state. The solid lines correspond to different Pade approximants of the two- (a) and one-loop (b) series obtained via the field- 
theoretical dimensional expansion [19114] , For d > 4 mean-field results become exact. The available numerical estimates obtained 
via Monte Carlo simulations [9120114] are also indicated with the corresponding errorbars, if provided. The actual agreement 
between analytical predictions and numerical results is quite good. 



force F ~ —5J-/5L acting on the confining walls is the 
sum of a bulk term -Fbuik and the Casimir force Fq which 
is actually the force acting on the confining walls of the 
film when they are "immersed" in the fiuctuating medium 
(so that the contribution i^buik acting on the two opposite 
sides of each wall cancels, see Fig.|4]). Note that this way of 




bulk 



Fig. 4. Fluctuation-induced forces acting on two parallel walls 
immersed in a fluctuating medium. The net force on each wall 
is given by the Casimir term Fc- 



determining the Casimir force Fc via the equilibrium free 
energy J- gives no information about the force at a given 
point in space. A different approach which, in principle, 
provides such an information is based on the stress tensor 
Tq^(x). Tq^(x) is a local functional of the order parame- 
ter field 'p(x), it is determined - in equilibrium - on the 
basis of 7i, and its expectation value on the equilibrium 



distribution function Pcq yields the force 
F = A:BTA(Tii(x))eq, 



(10) 



where T±± is a suitable component of the stress tensor Tap 
(see Ref. [25] for details). In equilibrium such an average 
is actually independent of the position x at which it is 
calculated; nonetheless it is tempting to interpret Ti^(x) 



as providing the local fluctuating pressure (per /cbT) at 
the point x in space. Summing up, the equilibrium value 
of Fc can be determined via both the free energy J- and 
the stress tensor Tap{x.). On the other hand, if one is in- 
terested in relaxation phenomena the definition via J-' is 
no longer viable. Before discussing this point let us illus- 
trate a simple experimentally relevant setting in which the 
relaxation of Fc might be observed. Indeed, consider the 
case in which the Casimir force, in competition with dis- 
persion forces, determines the equilibrium thickness L of 
a wetting layer of a binary ^He/'^He mixture at a given 
temperature T (see, e.g., Refs. [26127) ), as in Fig.jSjJa). 



vapor 



] — ^He-'^He fiquid" T 
/ / P// /////// / / 

(a) 




Fig. 5. (a) The equilibrium thickness L of a wetting layer of 
a '^He/*He mixture is affected by the Casimir force, as experi- 
mentally demonstrated in Ref. [55]. (b) After a sudden change 
of the temperature from the initial value Tin to Tfln, the Casimir 
force is expected to change in time so that the fllm, originally of 
thickness Lin, attains its final thickness Lfln after a relaxation. 



Measuring L as a function of T allows a determination 
of the equilibrium Casimir force. Upon changing suddenly 
the temperature (and therefore the correlation length ^ of 
the concentration fiuctuations in the mixture) the Casimir 
force Fc changes and the thickness of the layer relaxes to- 
wards its final value (as sketched in Fig. [5t^b)). Although 
experimental data for the evolution of L are not currently 
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Fig. 6. Universal scaling function ■3R{t) of the Casimir force Fcit) in d = 3 after a quench from the high-temperature phase 
to the critical point (see Eq. (|13|) ) in the Gaussian model with purely relaxational dynamics. The finite-size scaling variable 
associated with time is t = {t/TR){(,o/L)^ . The scaling function has been normalized by its asymptotic (equilibrium) value 
•&R{t oo) (see Eq. (|12|) '). (a) i?H(t) for symmetric boundary conditions (DD and NN), leading to an attractive Casimir force 
which is equal on both the confining walls, (b) i}R{t) for the case in which the walls impose difi'erent boundary conditions on the 
order parameter (DN or ND). In contrast to equilibrium, the dynamic value and the qualitative features of the time dependence 
of Fc depend on the boundary condition imposed by the wall at which Fc is measured. 



available, some aspects of the dynamics of fluctuation- 
induced forces are expected to be accessible in the near fu- 
ture. The evolution of L is coupled to that one of the order 
parameter field (and other relevant hydrodynamic quan- 
tities) which evolves in time under the effect of thermal 
fluctuations; in turn, due to universality, this dynamics 
can be studied, close to critical points, via simplified mod- 
els (typically in the form of Langevin evolution equations 
for the order parameter), as already mentioned in Sect. [5] 
(see Ref. [2] for review) . Instead of discussing the problem 
within the experimental setup described above, in which 
L changes in time, we assume that the fluctuating medium 
is confined between two parallel walls at fixed distance L 
while the force Fc{t) acting on them is measured (Fig. 2]). 
As in the case of static behaviour, the presence of the 
confining walls imposes boundary conditions on the time- 
dependent order parameter field ip{x,t), the dynamics of 
which is affected by the confinement. In order to determine 
the dynamics of the Casimir force one has to understand 
how the local order parameter translates into the local 
force which can be measured, say, at the walls. One natural 
connection between the two is provided by the local stress 
which, via the time-dependence of the order parameter, 
becomes time dependent: Tq^(x, t) = Taf3\ip(^x)i-^ip{x.t)- Ac- 
cordingly, it is natural to define the dynamic force acting 
on the walls W as [28] (compare with Eq. (fTO)) ) 

F{t) = kBTA{T^^{xeW,t)) noise (11) 

where the average (. . .)noisc is taken over the possible real- 
ization of the thermal noise. Heuristically, this amounts to 
the assumption that at each time there is an "energy cost" 
kBTA5L{T±±) associated with the displacement 5L of one 
of the confining wall, which is actually determined, as in 
equilibrium, by the order parameter field. Note that in 
thermal equilibrium this definition renders the static force. 
(See Ref. [22] for some proposed definitions of fluctuation- 
induced forces in different non-equilibrium system.) At 



variance with the equilibrium case, however, the local pres- 
sure F/A might actually depend on the point x at which 
it is measured. In what follows we assume translational 
invariance in the direction parallel to the walls, yielding 
a spatially constant F on each of them. Finally, in or- 
der to obtain the dynamic Casimir force Fc{t) one has to 
subtract from F{t) the corresponding bulk contribution 
Fhu\k{t) = limL~*oo F{t). To see this definition at work 
we consider the case in which Tin 3> Tc, corresponding 
to a correlation length so small that each of the walls 
cannot feel the effect of the confinement imposed by the 
other, resulting in Fc{t < 0) = 0. For the final tem- 
perature, instead, we assume Tfin — which eventu- 
ally leads to a non-vanishing, long-ranged Casimir force 
Fc{t > 0) 7^ 0. (Note that, due to the critical-point shift 
in the film geometry [30j , Tc is not the critical temperature 
of the film.) The analytic results illustrated below refer to 
the case in which the confining walls impose Dirichlet (D: 
ip{x e W,t) = 0) or Neumann (N: d±Lp{x e W,t) = 
where d± is the derivative in the direction normal to the 
wall) boundary conditions onto the order parameter field 
(fi. We assume that the effective Hamiltonian H which cap- 
tures the equilibrium properties of the system is quadratic 
in the order parameter (Gaussian model) and that the dy- 
namics is purely dissipative (Model A). (We point out, 
however, that the actual dynamic of fiuids and binary 
mixtures requires accounting for features which are not 
present in Model A [5]). As a consequence of the quadratic 
Hamiltonian, the relevant component of the stress tensor 
can be expressed, in terms of the order parameter field, as 
T±±ix,t) = {d±(p{x,t))^/2- (V||^(x,t))V2 where V|| is 
the gradient parallel to the wall. 



In the long-time limit the film attains equilibrium and 
the force takes different values depending on the combina- 
tion of boundary conditions imposed by the two confining 
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walls: 



Fc 
A 



A 



(„s) 



< 0, for DD, NN, 
> 0, for DN, ND, 



(12) 



where, in three dimensions d = 3, the amplitudes for 
symmetric and non-symmetric boundary conditions are 



A 



0.048 ... and Z\^"'^ = +0.036 . . ., respectively 
(see, e.g., Ref. [2S])- According to Eq. (fT2|l . the Casimir 
force leads to attraction between the confining walls if they 
impose the same boundary conditions on the fluctuating 
order parameter (NN, DD), whereas repulsion results from 
different boundary conditions (DN, ND). Note that in the 
latter case the equilibrium force Fc acting on the two 
walls is the same, independently of the fact that they im- 
pose different boundary conditions on the field; in this 
sense they cannot be distinguished by measuring Fe- 
in view of dynamical scaling, the relaxation of the force 
acting on the confining walls after the thermal quench is 
expected to depend on the temporal scaling variable t — 
{t/To){^o/L)^ via the dynamical scaling function i9_R(t), 



Fc{t) fesT, f t f^o 



(13) 



where z = 2 is the dynamical exponent of the model we 
are focussing on. At small and large times we recover the 
initial vanishing force and its final equilibrium value in 
Eq. (fT2|) . respectively, yielding i?_r(0) — and ??_R(t 



oo) 



A 



(s,ns) 



depending on the actual boundary con- 



ditions. In Fig. [6] we report the analytic predictions for 
'0R(t) (normalized to the corresponding equilibrium value 
'dji{t — > oo)) for the case of symmetric DD, NN (a) and 
non-symmetric DN, ND (b) boundary conditions. The ap- 
proach to the asymptotic value depends on the boundary 
conditions but whereas in (a) it is independent of the wall 
at which Fc is measured (due to the symmetry), this is 
no longer the case in (b) . Accordingly, the dynamics of Fc 
provides a way of determining which one of the two walls 
imposes D or N boundary condition in the DN (or ND) 
case, even though this is not possible by simply looking 
at the equilibrium force. Interestingly enough, the time 
dependence is not monotonic in the case of the D wall. 

The definition pT|) allows one to study the effect of an 
external perturbation on the Casimir force acting on the 
confining walls, as discussed in detail in Ref. within 
the present model and DD boundary conditions. In par- 
ticular, if the external field (conjugate to the order pa- 
rameter) is localized in time and space, as indicated by 
the black dot in the side view in Fig. [Tlja), the response 
starts to propagate in the film, yielding asymmetric ef- 
fects. At the very early stages, the response has practically 
not yet reached the confining walls so that the force acting 
on them is basically the equilibrium one corresponding to 
a vanishing external field. In course of time the pertur- 
bation induced by the field hits the confining walls and 
correspondigly the force exerted on them increases (re- 
ducing the equilibrium attraction), with different magni- 
tude at different points. Finally, because of the relaxation 



L/3 



(a) 



(b) 



Fig. 7. Effect of an external field on the Casimir force act- 
ing on the confining walls. The impulse field is applied at time 
t — and distance L/3 from the left wall (a). The response of 
the medium is such that the Casimir force locally increses com- 
pared to its (negative) equilibrium value, reaches a maximum 
and then, due to the relaxational character of the dynamics, 
decreases again to its equilibrium value. The spatial depen- 
dence of the total local force at time t = 0.096 for the critical 
Gaussian model with DD boundary conditions is sketched in 
(b), where the colorcode indicates increasing force from bottom 
to top, compared to the equilibrium value. The magnitude of 
the field-induced effect depends on the square of the strength 
of the applied perturbation (see Ref. [28] for details). 



character of the dynamics, the perturbation induced by 
the external field vanishes in the limit of long times and 
the effective force reaches again its equilibrium value. In 
Fig. [Tljb) we report a sketch of the distribution of the re- 
sulting stess on the confining wall at time i — 0.096 in 
the case in which the impulse has been applied at time 
i = in the position indicated in Fig. Elja). The system is 
assumed to be at the bulk critical temperature Tc- Due to 
the different distances from the two walls, the maximum 
force for a; = is attained earlier on the left wall com- 
pared to the right one. (We refer the reader to Ref. [28] 
for further details.) 

In summary, we have proposed a definition of the dy- 
namic Casimir force Fc via the local and time-dependent 
stress tensor (see Eq. (fTTjl). which allows one to study 
the equilibrium and non-equilibrium dynamics of Fc in a 
variety of different cases, dynamics and boundary condi- 
tions. The analysis of very simple models reveals already 
a quite rich behaviour. However, it is desirable to estab- 
lish a clearer connection between this definition of the dy- 
namic Casimir force and the force that can be measured 
directly in actual experiments and molecular dynamics 
simulations. 



4 Perspectives 

In the previous Sections we have reviewed some aspects 
of relaxation phenomena in bulk and confined critical sys- 
tems. Here we mention some issues which in our opinion 
deserve further investigation. 

It is well known that surfaces affects locally both the 
static and dynamic critical behaviour |24|31j and there- 
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fore we expect that also non-equilibrium properties such 
as the scaling behavior of two-time response and correla- 
tion functions after a quench is accordingly changed. In 
spite of some available Monte Carlo data and preliminary 
analyses |32lll| , a detailed and quantitative study within 
the field-theoretical approach is still lacking. The slow 
collective relaxation phenomena discussed in Sect. [2] are 
due to the divergence of the equilibration time upon ap- 
proaching an equilibrium critical point (characterized, e.g., 
by detailed balance). A qualitatively similar behaviour, 
however, has to be expected also in the cases in which 
the system undergoes a non-equilibrium phase transition 
of dynamical nature (e.g., reaction-diffusion systems, see 
Ref. [33j and references therein). Also in this case, field- 
theoretical methods provide valuable insight into the prob- 
lem [33] and it would be desirable to extend the analysis 
to other relevant non-equilibrium universality classes. A 
particularly intriguing question is how (and if) the non- 
equilibrium critical relaxation phenomena in the bulk chang( 
when quantum fluctuations come into play, e.g., upon ap- 
proaching a quantum phase transition point. 

The study of the dynamics of fluctuation- induced forces 
is still far from being satisfactory and complete. In order 
to allow a comparison with experimentally relevant set- 
tings (e.g., involving colloidal particles in suspension) it is 
particularly important and urgent to explore the effects of 
geometries and boundary conditions different from those 
discussed here and to consider dynamics which are more 
suitable for the description of binary fluid mixtures, on 
which most of the experiments are based. 

I am grateful to P. Calabrese, S. Dietrich, F. Krzakala, A. Ma- 
ciolek, R. Paul and G. Schehr for the stimulating collaborations 
which lead to some of the results presented here. 
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